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Abstract 

It is shown that the observation of the spectra of leading neutrons from proton beams 
can be a good probe of absorptive and migration effects. We quantify how these effects 
modify the Reggeized pion-exchange description of the measurements of leading neutrons 
at HERA. We are able to obtain a satisfactory description of all the features of these 
data. We also briefly discuss the corresponding data for leading baryons produced in 
hadron-hadron collisions. 

1 Introduction 

Leading neutrons, which are produced from proton beams, are of special interest since the 
production process is dominated by tt exchange, see Section 2. Such leading neutron spectra 
have been measured recently in photon-proton collisions at HERA These data supple- 

ment measurements of leading baryon spectra made in hadron-hadron collisions many years 
ago [0]-jn]. By observing leading neutrons we have essentially a tagged n beam. For example 
at HERA this 7r beam may participate in deep inelastic scattering (DIS) or photoproduction. 
In this way both the structure function of the pion F£ |H 12] and the form of the underlying 
Reggeized 7r trajectory [3] can be studied. 

There is another aspect of such processes which merits particular study. Soft rescattering 
effectively leads to the absorption of leading neutrons with Feynman x (which we denote xl) 
close to 1; the role of this correction was originally studied in Refs. |11[ IT2"] . At that time, 
these were essentially theoretical studies, which provided predictions for future measurements. 
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Figure 1: (a) The pion-exchange amplitude and (b) the corresponding dominant triple- Regge 
contribution to the cross section of the inclusive production of leading neutrons, ap — > Xn. 
The coupling of the Pomeron— n — n Regge trajectories is denoted by r P7Tn (t) in (JTJ). 



Now detailed leading neutron data have become available from the experiments at HERA. 
These experiments measure the ^*p — > Xn cross sections for producing leading neutrons with 
different xl and pt values, as functions of both the virtuality, Q 2 , of the photon, and the 
cm. energy, W, of the incoming jp system. The HERA data offer the opportunity of seeing 
how absorption chang function of all these kinematic observables. We also show the 

importance of migration, or change of the kinematic variables, of the leading neutron due 
to rescattering effects. The ways to estimate the effects of absorption and of migration are 
discussed in Sections 3 and 4 respectively. The description of the HERA data, jp — > Xn, is 
presented in Section 5. For completeness, in Section 6, we discuss the behaviour of leading 
baryons produced in hadron-hadron collisions, in particular in the processes pp — *> Xn and 
pn — > Xp. 

In exclusive reactions, the rescattering could produce new secondaries which has the effect 
of suppressing the rate. That is, an exclusive cross section is suppressed by absorptive effects. 
On the other hand, for an inclusive process the rescattering, will just change the energy and 
transverse momentum of the leading baryon, and will depopulate the region of large xl and 
small pt- From the viewpoint of an experimental trigger, for xl — > 1 this appears as absorption, 
but it is better to consider it as a migration into an enlarged phase space. In particular, such a 
migration can affect the form of the effective (that is measured) 7r trajectory, leading to a larger 
mean pt at lower xl- Of course, this may also modify the measurement of the pion structure 
function, F%, which is extracted from the HERA data. 



2 Leading neutrons from n exchange 

To begin, let us recall the triple- Regge formula which describes the inclusive processes ap — > Xn 
with a = 7 or p. For xl close to unity the inclusive leading neutron spectra is given by the 
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triple-Regge diagrams shown in Fig. ^ Thus we have 



flT V s o/ sm^vra^t)) 

where t = (p p — p n ) 2 ', and gp(0) is the coupling of the Pomeron to the photon, via a qq pair, 
or the proton, depending on whether a = 7 or p. The above form, with only the dominant 
pion-exchange contribution, has been used by many authors, see, for example, Refs. [21 EH EI| • 
Moreover, both the old hadron-hadron interaction data and the recent HERA results on the 
spectra of neutrons |3] clearly indicate the reggeization of pion exchange. In fact, detailed 
studies of hadronic interactions in the reggeized pion-exchange model fBl E] have shown that 
the formula can be used in a broad range of xl, with xl ^ 0.4 — 0.5. 

The contribution of reggeized pion exchange to the inclusive production of baryons, ap — ► 
XN with N = p, n, is of the form 



da{ap -> XN) _ Gl +pn (-t) 



F\t) <C(M 2 ) (1 - xl) 1 - 20 *®, (2) 



dxLdt 167T 2 (t — m 
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where a n (t) = a' n (t — m 2 ) is the pion trajectory with slope a4 ~ 1 GeV -2 , and Gv 0pp /47r = 
G 2 + /87T = 13.75 15 j. The invariant mass M of the produced system X is given by M 2 ~ 
s(l — Xi). The signature factor, and the difference between sin(7ra 7r (t)) and (t — m 2 ) factors 
in the denominators of (0) and (J2J), are absorbed in the effective vertex form factor F(t). This 
form factor is usually taken to be of the form 

F(t) = exp(6t/2), (3) 

where, from data at relatively low energies, we expect b ~ 4 GeV -2 , see, for example, Ref. [T4*j . 
However more recent analyses find b ~ [TT)| Ej . 

For the pion-exchange contributions to the spectra of both protons and neutrons, it is 
possible to take account, not only of the diagram of Fig. |2k, but also, of the diagram of Fig. |2fc>. 
Besides this we can include the contributions of p, 02, ... exchanges, and also of resonance decays, 
such the A-resonance contribution in the ttN amplitude, to the leading p and n spectra. The 
calculations of the spectra of neutrons of Refs. O EJ EH], show that the diagram of Fig. |2K 
dominates, and that the contribution of the diagram of Fig. |2b is about 20% at most. On 
the contrary, the leading proton spectrum is dominated by the triple-Pomeron diagram; the 
diagram of Fig. |2t gives only about 30% of the inclusive cross section for x p ~ 0.8 in the ISR 
energy range. 

One interesting application of the leading neutron spectra at HERA is the possibility of 
using (J2J to extract information on the 77? total cross section, 0"*°*, for both real and virtual 
photons. This, in turn, allows a measurement of the pion structure function F£(xb, Q 2 ) at very 
high energy or small Bjorken xb- The direct application of (J2J) to the HERA photoproduction 
data j2] leads to the result a^/a^ = 0.32 ± 0.03, which is about a factor 2 lower than the 



3 




(a) 



(b) 



Figure 2: Diagrams (a) and (b) show the reggeized 7r-exchange contributions to the leading 
ap —>■ XN and ap — > X(Nn) spectra respectively, with N = p,n. 



expected 1 ratio of 2/3. However, so far, we have neglected the absorptive corrections, or so- 
called rescattering effects, to (J2J). These are very important at high energies. They reduce the 
predicted cross section, as well as modifying both the energy and the Q 2 dependence of the 
inclusive spectra of leading neutrons. We discuss the form of the correction in the next Section, 
and present the leading neutron spectra at HERA in Section 5. 



3 Unitarity effects and gap survival 

The Born-type regge pole exchange diagram of Fig. ^ represents only the first approximation 
to the description of the leading neutron spectra at high energies. Multi-pomeron exchanges 
introduce modifications of the amplitudes, and lead to the restoration of unitarity jTHJ EI] . 

We use the space-time diagram of Fig. El to distinguish between two types of rescattering. 
First we have the eikonal scattering of the leading hadrons (either incoming or outgoing) indi- 
cated by exchange (a). The other possibility is the rescattering of the leading hadron on one 
of the intermediate partons in the central rapidity region, indicated by exchange (b). 

In terms of Feynman diagrams, an additional Pomeron exchange leads to a negative contri- 
bution; that is absorptive corrections diminish the size of the predicted cross section. From the 
physical point of view, the inelastic interaction, shown by either exchange (a) or exchange (b) 
in Fig. El produces new secondary particles which populate the rapidity gaps and carry away 
energy from the leading neutron. In other words a leading neutron, with a large energy fraction 
XL-, can only be observed in the small fraction of events which have no secondary inelastic inter- 
actions. The rescattering corrections therefore reduce the probability that the leading neutrons 
will be found in the large xl bins. 

1 Note that experimental estimates of ct*°* at high energies, based on the absorptive corrections to the 
amplitude for 7 + p — > n + ir~ + p and its interference with the 7 + p — > p° +p amplitude, are about §0^°* at the 
same energy ^7]. Thus we expect the additive quark model estimates to be reliable. 



4 



y or p 




Figure 3: The space-time diagram for the amplitude describing leading neutrons produced by 
the inclusive processes 7p — > Xn or pp — > Xn. The two types of multi-Pomeron absorptive 
corrections are indicated symbolically by the shaded areas (a) and (b). The corresponding 
corrections to the cross section are shown in Fig. E]and Fig. El respectively. 
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Figure 4: Symbolic diagrams of the eikonal absorptive corrections to the cross section for the 
inclusive process ap — *> Xn. The extra lines denoted by P, which surround the triple-Regge 
interaction, represent multi-Pomeron exchanges between the leading hadrons. 
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Figure 5: Symbolic diagrams for the "enhanced" absorptive corrections to the cross section for 
the inclusive process ap — > Xn, which become important at very high energies. The mass M 
of the produced system X is assumed to be sufficiently large for both the rapidity intervals y\ 
and yi to accommodate Pomeron exchange. The extra lines denoted by P, which are coupled 
directly to the ingoing p or outgoing n, represent multi-Pomeron exchanges. 

Technically the absorptive corrections are calculated in the following way. The most familiar 
are the eikonal-type corrections which are shown symbolically in Fig. 0J where only elastic 
rescatterings are taken into account. The corrections are well known: 



where px is the impact parameter, and dao is the contribution to the cross section for leading 
neutron production from the lowest order diagram of Fig. db). The opacity Q is the impact 
representation of the single-pomeron contribution to elastic scattering, which can be determined 
from studies of the total and elastic processes; for a recent review see, for example, Ref. [20]. The 
formalism can easily be extended to take into account diffractive excitations of the initial and 
final hadrons. It leads to a multi-channel version of (j3J), see, for example [2IJ HE1 H23 • Diagrams 
of this eikonal type have been considered in previous estimates of absorptive corrections to the 
spectra of leading baryons [HJH21- 

Note that since Q(s,pt) increases with energy, the suppression of the Born cross section dao 
also increases with energy. The damping factor, exp(— Q(s, pr)), is strongest for small values of 
the impact parameter p?. Thus the absorptive suppression is greater for the p, a% reggeons than 
for 7r-exchange, which, due to the l/(t — m^) behaviour at small t, has a broader distribution 
in p T . 

The unitarity corrections have an interesting effect on the t-dependence of the cross section 
driven by pion-exchange, (j2J. As the pion pole is approached, that is as t — > the corre- 
sponding impact parameter becomes very large, pr — > oo, and the gap survival damping factor 
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Figure 6: Multi-pomeron corrections to the reggeons in the triple-regge diagrams. They have 
no, or at most a weak, dependence on energy. The curved lines denoted by P represent multi- 
Pomeron exchanges. 



Thus we anticipate a more complicated t dependence than that given by eqs. (J2J) and (jHJ). First, 
at very small t, \t\ < 4m 2 , the survival factor S 2 rapidly decreases from S 2 = 1 at t = m 2 , 
finally reaching a value ~ exp(— Q(pt = 0)) at large \t\. However before then, the additional t- 
dependence coming from the pion vertex and trajectory reveals itself, namely fv(t) = exp(bt/2) 
with b ~ 1 GeV 2 . So, in (J2J), we have 

F(t) = f v (t)S(t). (6) 

Thus fitting the data using the single exponential form Q may confuse the situation. Note 
that in the discussion below we will use the average value of S* 2 

S 2 = <exp(-fi(p T ))>. (7) 

Another class of multi-pomeron diagrams corresponds to the modification of the upper 
reggeon in the triple-pomeron diagram (Fig. Ek)> or of the lower reggeons (Fig. IHt 1 ) , or of both 
(Fig.Efc). The diagrams of Fig. [H^l modify only the M 2 dependence of the cross sections, and for 
the 7r-exchange model can be absorbed in the behaviour of a^(M 2 ). The diagrams of Fig. |Ub 
modify the xl dependence, and for the 7r case can be by modifying the 7r-trajectory (essentially 
the slope a' n ) and the form factor. It is important, that to a good approximation, they do not 
depend on the energy 2 y/s. 

The multi-pomeron diagrams of Fig. are only expected to become relevant at sufficiently 
high energies. That is when the energy is large enough for the mass M of the produced system 
X to itself be large enough for both the rapidity intervals y\ and y 2 = Um — yi to be sufficiently 
large to accommodate pomeron exchange 3 ; that is y$ > yo with yo — 2 — 3. Recall that 

2 The unfactorizable diagram Fig. [§p can, in principle, have some weak s dependence. 

3 Diagrams of this type were mentioned in our paper on dijet production at the Tevatron but for this 
hard diffractive process they were small even at Tevatron energies. 
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Figure 7: A cut leading to an s-channel discontinuity of the triple-pomeron diagram for the 
elastic ap scattering amplitude, 

y M = ln(M 2 /M 2 ) with M 2 = 1 GeV 2 ; so for a typical value of M 2 ~ 3000 GeV 2 we have 
Vm = Vi + Jj2 ~ 8. Note that rapidity intervals yi < yo are included in the contributions of the 
diagrams of Fig. H]and Fig. Efb). 

Let us study in more detail the diagram shown in Fig. Et, which describes the process 
ap — > Xn. This diagram is one of the s-channel discontinuities of the triple-pomeron diagram 
for the elastic ap scattering amplitude 4 , see Fig. |7J Indeed the diagram is related to the triple- 
pomeron diagram for the elastic amplitude by the AGK cutting rules P3j. The corresponding 
coefficient is 4. So it is straightforward to obtain the correction to the ap — * Xn cross section 
from diagram of Fig. provided we know the cross section of Fig. [7| in the triple- Pomeron 
region, xl —> 1. By fitting this cross section to the triple-Pomeron formula we can extract the 
strength of the vertex rppp [22] • If we take the value of the effective vertex rppp determined in 
the analysis of the old hadron data, then we find the correction due to the enhanced absorptive 
diagrams to be about 15%. 

So there is a relatively small contribution coming from the enhanced graphs of Fig. 
that is from the rescattering of intermediate partons, as indicated by the shaded region (b) of 
Fig. El This is consistent with the fact that within the HERA range (W = 40 — 270 GeV) no 
energy dependence is observed in the leading neutron yields; see, for example, Figs. 6 of [4 5 , 
and Tables 14, 18 and Figs. 11, 12 of [2] which show, for fixed Q 2 , the same probability 6 to 
observe a leading neutron for values of x Bj which decrease by more than an order of magnitude 
corresponding to an increase of the photon laboratory energy by more than a factor of 10. This 
flat experimental behaviour of the probability indicates that actually the absorptive corrections 
caused by the rescattering of intermediate partons are much smaller than might be expected 

4 The summation of more complicated enhanced diagrams was discussed in [52], and more recently in |23| . 
5 See also Figs. 2 and 4 of |Sj for leading proton data, where the effect of enhanced diagrams should be the 
same. 

6 That is the same probability, S 2 , to observe the rapidity gap associated with pion-exchange. 
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from leading-order perturbative QCD calculations, similar to those presented in Ref. (see 
also the discussion in (23 |2E| ) ■ 



4 The inclusion of migration 

In order to compute the spectra of leading neutrons, we must consider the effects of migration, 
as well as of absorption. We assume that each rescattering spreads out the Pt,xl spectra of 
the leading baryon according to distribution 7 

(1 - a)xl a b m exp(-p|6 m ), (8) 



dp T dx l 

where we take the slope b m = 6 GeV -2 , and the intercept a ~ — \ corresponding to secondary 
Reggeon exchange in the Kancheli-Mueller approach. Thus, after two rescatterings 

^ (2) 1.2 f d 2< lT , 2i. \ ( f- - \2z. \ f 1 dx±dx 2 M N2 :/ \ 

b m / exp(-g T 6 m ) exp(-Q» r - q T ) b m ) I —{I -a) 5{x L -xiX 2 ) = 



dp^dXL m J 7T ^ V ^ 1 " l> L y ^ ± ' " lJ Jq (XiX2) a 

exp(-p 2 T b m /2) . (9) 



a) 2 ln(l, './•,.) ., 



2(x L ) a 

Correspondingly, after k rescatterings the distribution becomes 



dNW b m {l-af \n k -\l/x L ) 2 

- exp(-p T b m /k) . (10) 



dp\dx L k(k — 1)! (x L ) a 

The probability of rescattering, or, equivalently, the mean number v of rescatterings, at fixed 
impact parameter p T is given by the opacity v = ft 

%) = exp(-p|/4 J B r ) . (11) 

Here we focus attention on migration effects in jp Xn, so px = P^n- The cross section a 
is, however, not the total cross section, since the process proceeds in two stages. First the 
photon fluctuates into a qq pair, which may be considered as a sum of vector mesons, and then 
the vector mesons rescatter on a nucleon. We therefore assume a is the np total cross section 
having in mind the additive quark model and/or p meson dominance. B r is the slope of the 
rescattering amplitude (which we take to be B r = 5 GeV -2 ). 

Note that the impact parameter px in (jllj) . which controls the absorption and migration 
effects, is not equal to the parameter p nN which describes the space structure of 'pure' one-pion 



7 Strictly speaking, besides the migration in pt,xl space, there may be charge- exchange (that is neutron to 
proton transitions). We neglect this effect since the goal of the present section is not a precise quantitative 
description, but rather a qualitative evaluation of the role and possible size of the 'migration' phenomena. For 
the same reason we use a simplified form, of the leading baryon distribution. 
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exchange in leading neutron production. The subscript it denotes the 'transverse' position of the 
7i — 7i— Pomeron vertex in Fig. ^b). That is p n ^ is Fourier conjugate of the neutron transverse 
momentum in the bare (Born) amplitude leading to the cross section in (J2J). However, due to 
the rather large values of p n N, we may expect that the absorptive and migration effects will not 
be too strong — since the transverse distance between the leading baryon and the incoming 
photon (or gg-pair), 

PT = PyN = \PnN + Pyn\, (12) 

is, in turn, relatively large and the mean number of rescatterings v(pt) — ^{pt) at the periph- 
ery of the interaction (described by the amplitude fl) is rather small. Here p 77r is the impact 
parameter for the amplitude describing the interaction of the incoming photon with the 'effec- 
tive' pion (exchanged in the t-channel); recall that the subscript tt denotes the position of the 
7r — 7i— Pomeron vertex in Fig. 0] 

To account for the fact that pt and ar e not the same, for each value of p n N, we calculate 
the probability 

w(pt,Pttn) = N(p 7rN ) J d 2 p yn ti^Pyx) 5(p T - |pV + Pitn\) , (13) 

where fi^p^) is the amplitude in impact parameter space for the photon-pion interaction. It 
is of the form (jllj) with the same slope B e g = B r = 5 GeV~ 2 . Since the expression (J2J) already 
includes the probability of the 'effective'-pion-7 interaction the probabilities w(pt,P-kn) are 
normalized to one; that is the normalization N^p^n) is fixed by the condition 8 



dpr w(p T ,p n 



N 



(14) 



for every value of p T ^. 

Thus, to calculate the inclusive cross section for leading neutron production we have to 
integrate over the impact parameter px 



da 



dx^dp, 



d 2 p T 



exp(ip T ■ Pkn] 
2tt 



d 2 pnN w(p T ,p 7rA r)J'(p r ,p 7rA r,X L ) 



(15) 



where 



F(Pt,PitN,xl) = v{pt)c 



-v(pr) 



F [ °\pnN,X L )+ / dx' L ^°\ PwN ,x' L ) ^\p nN ,X L /x' L ) 



J-L 



(16) 



8 Strictly speaking we should work with amplitudes and not with cross sections. Moreover the impact pa- 
rameter pt in the amplitude A may not be equal to the parameter p T in the conjugated amplitude A* . We find 
Pt — Pt only after integration over the neutron transverse momentun pt- However to simplify the computa- 
tions, here we use the semiclassical approximation, which assumes pt = Pt, using the Fourier transform just 
for the cross section. This approach, which is usually used in Monte Carlo simulations, may be justified by the 
fact that the effect of migration becomes important only at relatively large pr, see Fig. |H1 So after we integrate 
up to pt = 0.5 — 1 GeV we find that the equality pt — p T holds to rather good accuracy. Nevertheless, the fact 
that we need to correct the normalization in i|13ll4|) is just the result of this semiclassical approximation. 
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The factor z/(pr)exp(— v(pr)) accounts for the probability of the first interaction and the ab- 
sorption at a given impact parameter p?, where v(pr) = &{pr) is the mean number of inter- 
actions. The first term in the square brackets, J-(°\ is the original cross section (J2J) in the p w ^ 
representation, that is 

d 2 PT -> * d<T<® 



nN) ~ — 

r 



dxidp, 



The second term function J 7 ^ accounts for the migration. When we sum over k = 1,2,... 
rescatterings, the spectra in p wn representation takes the form 



dpipdxL ^ dpipdxL J 2ir ^ [k — 1)! (xlY 

/j2 
Y~ exp(ip T • p) - a)e-" 2/4b ™ exp (v{\ - a) e ~ p2/4bm ln(l/x L )) /{x L ) a . (17) 

Thus the function T^ r ' in (jTSj) is 

^ r \p,x L ) = u(l-a)e- p2/4bm exp (u(l - a) e ~ p2/4bm ln(l/x L )) /(x L ) a . (18) 

The results shown below correspond to a = 0, but the spectra obtained for a = | are very- 
similar. 



5 Predictions for leading neutrons at HERA 

To predict the absorptive corrections to the pion-exchange formula in which we set 6 = 0, 
we use (HJ), together with (JTTjl in which we take a — 1.3 x 31 mb. The value o" to t( 7r p) = 31 mb 
agrees with the Donnachie-Landshoff parametrization [2H] and with the cross section evaluated 
by ZEUS via the 7r-proton absorption in 7p — > (7r + 7r~) p process The factor 1.3 takes 

account of the diffractive excitations of the initial and final hadrons [30j. 

The resulting predictions for the xl spectra of leading neutrons produced at HERA are 
shown in Fig. |SJ In the first plot the ZEUS acceptance cut, 9 n < 0.8 mrad, has been im- 
posed. In the next three plots we show the effect of imposing different pt cuts, rather than 
the 9 n cut, on the leading neutron. Each plot shows three curves. The dotted curve is the 
Reggeised 7r-exchange prediction calculated from (J2J), whereas the continuous and dashed curves 
are computed from (|15J1 and (jlOj) with and without migration (JF^) included. The data points 
correspond to the leading neutron spectra measured by ZEUS 0. 

As anticipated, the migration practically does not affect the region of large xl > 0.8, where 
it is essentially enough to account for absorption only. However at smaller xl, the role of 
migration becomes quite important. First of all, it compensates the absorption (which occurs 
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Figure 8: The predictions for the xl spectra of leading neutrons corresponding to the ZEUS 
kinematics. The dotted, dashed and continuous curves are respectively the results assuming 
first only Reggeised 7r exchange, then including absorptive effects, and finally allowing for 
migration. The different plots show the effects of imposing different pt cuts on the leading 
neutron. The data points are the spectra obtained from the measurements of Ref. [2|. Note 
that the curves plotted here account only for the pion-exchange contribution, which is not the 
dominant source of leading neutrons for xl ^ 0.4. 
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for xl close to 1), and provides the conservation of baryon charge. Next, it changes the 
behaviour of the xl and pr distributions of the leading neutrons; and modifies the 'effective' 
slope of the reggeized pion trajectory which was extracted from the data using the simplified 
formula (0). 

From Fig. |H1 we see that the predicted leading neutron spectrum is in satisfactory agreement 
with that measured by ZEUS. The absorptive corrections reduce the cross section given simply 
by Reggeised pion exchange by a factor of just less than 0.5 so that the predictions are in 
agreement with the data at large xl- The reduction by 0.4-0.5 is not inconsistent with an 
earlier calculation [2] which gave a rapidity gap survival factor S 2 of 0.34 for the resolved part 
of the photon wave function. 

Moreover the data on the production of a leading neutron, together with either a pair of 
high Et jets jS2] or charm (D*) [33], are consistent with a much larger value of S 2 , as is to 
be expected when we select x 7 — > 1 events which sample a point-like photon, which directly 
interacts with high Et dijets. In particular, the probability to observe a leading neutron 
increases by up to a factor of 2.5 (that is S 2 increases up to about 0.9) when the momentum 
fraction carried by the dijet grows from to 1, see Fig. 8c of Ref. [32] ■ In Ref. [31] the HI dijet 
data were described, within NLO QCD [32], by pure pion-exchange without the inclusion of 
the survival factor, that is with S 2 — 1. However the t-dependence was not studied. Moreover 
a larger t-slope parameter, R = 0.93 GeV^ 1 , was used in [33] than the value R = 0.5 GeV -1 
used in |3~5] . Thus the b slope in [33] is 3.5 times larger. In this way the absorptive corrections 
were mimicked in [33], see the discussion leading to eq. (jUJ). 

The comparison with the ZEUS data in Fig.[S]shows evidence of migration for 0.4 < xl < 0.6. 
The points at lower x^ indicate that there may be additional contributions 9 caused by the re- 
combination of an initial valence quark with a sea quark in the neutron with xl ~ 0.3, as well 
as by baryon charge transfer to the central region 10 , whose description is beyond the present 
analysis. 

The leading neutron measurements at HERA also allow us to explore the Q 2 dependence of 
the problem. Unlike the naive Vector Dominance Model (VDM), where the effect of absorption 
should disappear as l/Q 2 , in QCD (or in the Generalized VDM) these effects are expected 
to decrease as 1/lnQ 2 . It was shown in Ref. [2] that the Q 2 behaviour observed at HERA is 
consistent with predictions [TT1 [T2] based on the QCD dipole approach. 

6 Leading baryon spectra in hadron-hadron collisions 

Some years before the HERA data became available, leading nucleon spectra had been observed 
in hadron-hadron collisions in a variety of experiments, albeit at lower energies. Generally, 

9 Note that in our computation we have already used b = 0. Thus we have no possibility to enlarge the 
neutron yield for xl ~ 0.5 by diminishing the value that we take for the slope b. 

10 The ratio p/p > 1 measured in the central region at RHIC and HERA indicates the presence of such an 
effect. 
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these, and related, data are well described by the triple- Reggeon formalism. At first sight, it 
appears strange that no gap survival factor suppression, S 2 , was necessary — particularly in 
view of our discussion of the HERA data in the previous section. However, this factor is effec- 
tively included in the normalization of the triple-Reggeon vertices [23] which were determined 
by fitting the data. The only problem is the 7T7r-Pomeron vertex, which must be consistent 
with the known result for the on-shell pion as t — > m 2 . The on-shell vertex is known from 
an independent analysis of np total and differential cross section data. Now, it is crucial to 
note that the 7T7r-Pomeron vertex depends strongly on the value of t — m 2 , leading to an ef- 
fective factor F(t) = 1 at t — m 2 . In the physical region of negative t this factor is expected 
to be S 2 ~ 0.4 in the CERN ISR energy range jSHl E3 ED]- Indeed, the indications for such 
absorptive corrections to 7r-exchange were observed in data for pp —>■ XA ++ production 8, 36J. 
Spectra of leading baryons in the Reggeized pion-exchange model have been described in papers 
[T31 EH1 IIHI [HH EH]- In papers [I3J EE1 EH] a phenomenological form-factor, which to some 
extent mimic absorptive effects, has been used. 

The situation with inclusive nucleon-charge exchange in pp (and pd) collisions, that is in the 
spectra of neutrons in pp-collisions and protons in pn-collisions, is not clear. The largest set of 
the highest energy data was obtained at the CERN ISR, Refs. || and [7] . These ISR spectra of 
neutrons have a rather different form. It was emphasized in Ref. |HB] that the data of [E] are 
strongly different 11 from predictions of the pion-exchange model, and lead to a violation of the 
energy-momentum sum rule, while the data of [7] are in a reasonable agreement with theoretical 
predictions [39J. As mentioned above, the models of Refs. [311133 are based on Reggeized pion- 
exchange, with the gap survival factor S 2 represented by a form factor which, however, gave 
only a rather small effect (~ 10%). The equivalent process pn — > Xp was also measured [2] 
over range of energies, which partly overlap with ISR energy range and are consistent with the 
model of Ref. [3H] ■ On the other hand there are data on the process pn — > Xp at relatively low 
energies (11.6 GeV), which are lower than predictions of the pion-exchange model without the 
S 2 factor. 

Thus the final picture for hadron-hadron processes is not clear. Some experiments are in 
agreement with the theoretical expectations of appreciable absorptive corrections to 7r-exchange, 
while others appear to be described by simplified 7r-exchange. 

7 Summary 

We have described all the main features of the leading neutrons observed at HERA, both in 
photoproduction and as a function of Q 2 . For photoproduction we may refer to Fig. [HI We see 
that Reggeized pion-exchange on its own is insufficient to describe the production of leading 
neutrons. However the inclusion of absorptive corrections are found to reduce the 7r-exchange 

11 Note that in the CERN ISR paper [Sj the cross section plotted in Fig. 7 is about 3 times smaller than that 
given in the Tables and in Fig. 6. 
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prediction by a factor of just less than 0.5, and to bring the theoretical expectation in line 
with the data. Thus, after accounting for the rapidity gap survival factor S 2 , such data can be 
used to measure the 77r cross section and the pion structure function. Moreover the value of 
the gap survival factor S 2 may be monitored by comparing the spectra of leading neutrons in 
dijet production at low x 1 with those for x n — ► 1, where the cross section is dominated by the 
direct, point-like photon. Next, the dependence of the survival factor S 2 (which specifies the 
probability to observe a leading neutron) on the photon energy may be used to experimentally 
probe the contribution of the enhanced absorptive corrections of Fig. [SJ The data appear to be 
flat in energy indicating only a small contribution from enhanced diagrams. This is relevant to 
the proposed Higgs searches in exclusive diffractive production at the LHC Thus we conclude 
that the additional suppression of the exclusive Higgs cross section caused by the enhanced 
diagrams cannot be as large as mentioned in Ref. |2E!> see [2H] for a detailed discussion. 

We also found that the neutron spectrum shows some evidence of migration effects, see 
Fig. |SJ Including migration slightly reduces the value extracted for the slope of the pion Regge 
trajectory; the reduction is about 0.05 — 0.1 GeV~ 2 depending on the interval of xl used to 
measure a' w . 

Finally, in Section 6, we briefly reviewed the available experimental data for leading baryons 
produced in hadron-hadron interactions, which were obtained about 30 years ago. Here the 
experimental situation is confusing. Some experiments are in line with theoretical expectations, 
while others are not. 
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